function [NewIRF,SOLUTION]=TempFun_NewIRF(PP,SS,MODEL,EquJac,StdVec)

EquJac.NK       =   Equ_NK(PP,SS,1);
EquJac.Portfolio=   Equ_Portfolio(PP,SS,MODEL,1,[],[],0);
fir_ord         =   JacAssemble(MODEL,EquJac);
[gx,hx,gxhx_ExitFlag,gxhx_Info]=gx_hx(fir_ord,1);
if gxhx_ExitFlag~=1
    warning('Perturbation Solution is problematic');
    NewIRF      =   [];
    SOLUTION    =   [];
    return
else
    solution    =   struct('gx',gx,'hx',hx);
    IRF         =   IRF_1order(MODEL,solution,StdVec,[],4*10);
    NewIRF      =   TempFun_AugIRF(PP,SS,MODEL,EquJac,IRF);
    SOLUTION    =   struct('PP',PP,'SS',SS,'MODEL',MODEL,'EquJac',EquJac,...
                           'fir_ord',fir_ord,'gxhx',solution,'IRF',NewIRF);
end